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Abstract. Possibilities and difficulties of applying the theory of magnetic field generation by 
, convection flows in rotating spherical fluid shells to the Giant Planets are outlined. Recent 
' ^ ' progress in the understanding of the distribution of electrical conductivity in the Giant Planets 
suggests that the dynamo process occurs predominantly in regions of semiconductivity. In con- 
^ , trast to the geodynamo the magnetic field generation in the Giant Planets is thus characterized 
by strong radial conductivity variations. The importance of the constraint on the Ohmic dissi- 
pation provided by the planetary luminosity is emphasized. Planetary dynamos are likely to be 
of an oscillatory type, although these oscillations may not be evident from the exterior of the 
■ planets. 
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> ■ 1. Introduction 

. Even before the Pioneer and Voyager missions evidence for a global Jovian mag- 
netic field had been obtained through the observations of decametric radio waves by 



Burke & Franklin (1955) The subsequent discoveries of the magnetic fields of Saturn, 
' Uranus and Neptune have been highlights of the exploration of the solar system by space 
probes. While Jupiter and Saturn possess magnetic fields which are dominated by a dipole 
part that is nearly aligned with the axis of planetary rotation, Uranus and Neptune ex- 
hibit magnetic fields that do not show such an alignment and are also characterized by a 
strong quadrupolar component. This latter property has been interpreted in earlier pa- 

o 



pers (Ruzmaikin fc Starchenko, 1991[ [Stanley fc Bloxham, 2004 1 as a result of the thin 



' shell to which the dynamo process is assumed to be confined. This explanation is not 



, entirely convincing, however, since Liu et al. (2006)' argue convincingly that the dynamo 



^ process in Jupiter and Saturn is also confined to thin fluid shells. As will be pointed out 
• ^ in section 5 magnetic fields without a dominant axis-aligned dipole are a rather common 
/\ feature of numerical simulations of planetary dynamos. 

^ . In the following we shall first briefly introduce the mathematical background for numer- 
" " ' ical dynamo simulations. At the present stage it appears to be appropriate to introduce 
only a minimum of physical parameters in order to obtain an overview of the parameter 
dependence of planetary dynamos. Hence we shall not even consider the compressibility 
of the fluid in the Giant Planets and restrict the analysis to the Boussinesq approxima- 



tion. An alternative approach has been used by Evonuk & Glatzmaier (2006) (see also 
article by Glatzmaier in these Proceedings) where a fairly accurate representation of the 
compressibility of Jupiter has been attempted. 
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2. Mathematical Description of Convection Driven Spherical 
Dynamos 

For the equations describing convection driven dynamos in the frame rotating with the 
angular velocity we use a standard formulation which has also been used for a dynamo 
benchmark (Christensen et al., 200ip . A more general static state will be assumed with 
the temperature distribution Ts = Tq — f3(fr^/2 + ATrir^^{l — r])^^ where 77 denotes 
the ratio of inner to outer radius of the spherical shell and d is its thickness. AT is the 
temperature difference between the boundaries in the special case /3 = 0. In the case 
AT = the static temperature distribution Ts corresponds to that of a homogeneously 
heated sphere with the heat source density proportional to the parameter /3. The gravity 
field is given by g = —jdr where r is the position vector with respect to the center of 
the sphere and r is its length measured in units of d. 

In addition to d, the time (P' jv, the temperature j^ad^ and the magnetic flux density 
vi^liQ)^!'^ I d are used as scales for the dimensionless description of the problem where v 
denotes the kinematic viscosity of the fluid, k its thermal diffusivity, q its density and 
[I its magnetic permeability. The Boussinesq approximation is used in that q is assumed 
to be constant except in the gravity term where its temperature dependence given by 
a = —{dg/dT)/g = const, is taken into account. The dimensionless equations of motion, 
the heat equation for the deviation O of the temperature field from the static distribution 
and the equation of magnetic induction thus assume the form 

V^v + B • VB + re - Vtt = p-^{dtv + v • Vv) + rk x v (2.1) 

V-v = (2.2) 
V^e + [Ri + ReVr^^il - Tjy^] r • V = P{dt + v • V)e (2.3) 

Prnidt + V • V)B + V X (Ay X B) = P„B • Vv (2.4) 

where k is the unit vector in the direction of the axis of rotation and where Vtt includes 
all terms that can be written as gradients. The Rayleigh numbers Ri and Re, the Coriolis 
parameter r, the Prandtl number P and the magnetic Prandtl number are defined 

. = ^. P=-. P.,.^f (2.5, 

where Aq is a typical value of the magnetic diffusivity A which we allow to vary as a 
function of the distance r from the center in contrast to the other material properties. 
Such a variation reflects the often significant variation of the electrical conductivity a 
with radius according to the relationship A = cr^^/i^^. While P = 1 is often assumed 
with the argument that all effective diffusivities are equal in turbulent media we keep P 
as a parameter. 

Since the velocity field v as well as the magnetic flux density B are solenoidal vector 
fields, the general representation in terms of poloidal and toroidal components can be 
used, 

V = V X (V$ X r) + X r , B = V x (V/i x r) + Vg x r . (2.6) 

By multiplying the (curl)^ and the curl of equation (6a) by r we can obtain two equations 
for $ and 5* which will not be given here explicitly (see, for example, Simi tev fc Busse, 2005[ ). 
The equations for h and g are obtained through the multiplication of equation (2.4) and 
of its curl by r, 

[2(1 - a){r - n) + a]\/^L2h = P^idM - r • V x (v x B)] (2.7) 
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[2{l-a){r-ri)+a]\7^L2g + 2{l-a)r-^dr{rg) F„[ati25-r • V x (V x (v x B))] (2.8) 
where the L2 is defined by 

In (2.7) and (2.8) a linear dependence A — Ao[2(l — a)(r — r,) + a] has been introduced 
such that A = Aq at the middle of the layer, r = + 0.5. 

Either rigid boundaries with fixed temperatures as in the benchmark case ( [Christensen et al., 2001] ) , 

$ = 5r(r$) = * = e = at r^ri = r]/{l-r]) and at r = r,, = (1 - r?)"\ (2.9) 

or stress-free boundaries with fixed temperatures, 

$ = 92^$ = a^(^'/r) = e = at r = r, and at r = r,,, (2.10) 

are frequently used. The latter conditions allow to cover numerically a larger region of the 
parameter space since the thin Ekman layers at the boundaries are nearly absent. The 
radius ratio 77 = 0.4 is often used in the simulations since it provides a good compromise 
for the study of both, the regions inside and outside the tangent cylinder. The latter is 
the virtual cylindrical surface touching the inner spherical boundary at its equator. For 
the magnetic field it is convenient to employ electrically insulating boundaries such that 
the toroidal component of the field vanishes there, while the poloidal function h must be 
matched to the function /i'"^^ which describes the potential fields outside the fiuid shell, 

g = h-h^''> ^drih-h^''^)^0 a.tr = r, and at r = r,,. (2.11) 

Alternatively an infinitely conducting inner boundary can be assumed, 

drg = h = at r = ri, (2.12) 

but computations for the case of an inner boundary with no-slip conditions and an elec- 
trical conductivity equal to that of the fiuid arc also often done. The numerical integra- 
tion of equations (2.3), (2.7), (2.8), (2.9) and (2.10) together with boundary conditions 



proceeds most often with the pseudo-spectral method as described by Glatzmaier (1984) 
and Tilgner "(1999)1 which is based on an expansion of all dependent variables in spherical 



harmonics for the 9, ^-dependences, i.e. 

$ = ^V;'"(r,t)P;"(cos6')exp{im0} (2.13) 

and analogous expressions for the other variables, 5*, 6, ft. and g. P™ denotes the asso- 
ciated Legendre functions. For the r-dependence expansions in Chebychev polynomials 
are used. 

It should be emphasized that the static state v = B = G = represents a solution 
of equations (6) for all values of the Rayleigh numbers Ri and R^, but this solution is 
unstable except for low or negative values of the latter parameters. Similarly, there exist 
solutions with B = 0, but v 7^ 0, 7^ 0, for sufficiently large values of either Ri or R^ 
or both, but, again, these solutions are unstable for sufficiently large values of Pm with 
respect to disturbances with B 7^ 0. Dynamo solutions, as all solutions with B 7^ and 
|B| cx for r — > cxd are called, are thus removed by at least two bifurcations from the 
basic static solution of the problem. 



3. Convection in rotating spherical shells 

Past research on convection driven dynamos in rotating, self-gravitating spherical shells 
has shown that the properties of convection in the absence of a magnetic field do change 
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quantitatively, but not qualitatively after the Lorentz force enters the force balance. It is 
thus important to understand the properties of non-magnetic convection. A rough idea 
of the dependence of the critical Rayleigh number Rc for the onset of convection on the 
parameters of the problem can be gained from the expressions derived from the annulus 
model (see recent review of |Busse, 2002| 

^c = S ( T-T^ ) ' itanern)ir~h-i, (3.1a) 



+ P 

Pt 



^^pj (r™tan(?„0^2-^ ujc^\jp—^j 2"^ (tan" &„./r„0^ (3.16) 

where refers to the mean radius of the fluid shell, = (r^ + 7'o)/2, and 9m to the 
corresponding colatitude, 9m = arcsin (rm(l — 77)). The azimuthal wavenumber of the 
preferred mode is denoted by rric and the corresponding angular velocity of the drift of 
the convection columns in the prograde direction is given by uJc/nic- 




Figure 1. Banana cells in a thin rotating spherical fluid shell cooled from within. Convection 
driven by centrifugal buoyancy is made visible by a suspension of tiny flakes which become 
aligned with the shear. Because the inner sphere is cooled, the outer heated, the experiment 
models the planetary situation where gr avity and temperature gr adient are reversed since only 
the product enters force balance. (After [Busse fc Carrigan, 1976[ ) 

While expressions (3.1) correspond to motions in the form of columns aligned with 
the axis of rotation, convection in the form of "banana cells" is realized in less rapidly 
rotating or thinner spherical shells as experimentally visualized in flgure 1. An ana- 
lytical theory of the "banana cells" which includes the differential rotation generated 



by their Reynolds stresses was first derived by Busse (1970) At low Prandtl numbers, 
namely for P < 10/^/T inertial convection in the form of equatorially attached cells 
becomes prevalent (jArdes et al. 19971 [Simitev fc Busse, 2003[ ) . Its name reflects the fact 



that it can be described as a small modification of certain inertial waves ( Zhang, 1994 



Busse & Simitev, 2004 1 A third form of convection is realized in the polar regions of the 
shell which are defined as the two fluid domains inside the tangent cylinder. Since grav- 
ity and rotation vectors are nearly parallel in these regions (unless rj = ri jvo assumes a 
value close to unity) convection resembles the kind realized in a horizontal layer heated 
from below and rotating about a vertical axis. Because the Coriolis force can not largely 
be balanced by the pressure gradient in this case, the onset of convection is delayed to 
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Figure 2. Convection in rotating spherical fluid shells for r = 5 • 10^, R — 2.7 • Rc, in the cases 
P = 0.5 (upper row, R — 1.5 • 10^) and P — 20 (lower row, R — 4, ■ 10^). Lines of constant Ur 
in the middle spherical surface, r = + 0.5, are shown on the right. The middle plots show 
streamlines, rd^/dip = const., in the equatorial plane. The plots on the right indicate lines of 
constant mean azimuthal velocity u<^ in the left halves and isotherms of in the right halves. 



higher values of R where convection outside the tangent cylinder has reached already 
high amplitudes. 

More important than its influence on the onset of convection according to relationships 
(3.1) is the effect of the Prandtl number on convection at finite amplitudes. Typical 
features of low and high Prandtl number convection are illustrated in figure 2. The 
columnar nature of convection does not vary much with P as is evident from the two 
plots on the left side of the figure. At Prandtl numbers of the order unity or less, - 
but not in the case of inertial convection-, the convection columns tend to spiral away 
from the axis and thereby create a Reynolds stress which drives a strong geostrophic 
differential rotation as shown on the right side of the figure. This differential rotation in 
turn promotes the spiral tilt and a feedback loop is thus created. At high values of P the 
Reynolds stress becomes negligible and no significant mean tilt of the convection columns 
is apparent in the middle plot of the lower row. In this case the differential rotation is 
generated as a thermal wind caused by the latitudinal gradient of the axisymmetric 
component of 0. 

The properties of convection are refiected in the averages of the kinetic energy densities 
of the various components of the velocity field which are defined by 

= ^(1 V X (V^ X r) p), = 1(1 V* X r (3.2a) 

4 = ^(1 V X (V$ X r) ^( = 1(1 V* X r (3.26) 

where the angular brackets indicate the average over the fluid shell and $ refers to the 
azimuthally averaged component of while $ is defined by <I> = $ — $. Analogous 
definitions hold for the magnetic energy densities where E, $ etc. are replaced by M, h 
etc., In figure 3 energy densities have been plotted for convection with and without 
magnetic fields. 
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Figure 3. Kinetic (thin symbols) and magnetic (thick symbols) energy densities of convection 
driven dynamos as function of P for t = 3 • 10*, R = 3.5 • 10®, t] = 0.4, a — 1 (implying A = Ao) 
and values of Pm as indicated in the boxes. The components Xp, Xt, XpXt (where X — E or 
M) are represented by circles, squares, plus-signs and crosses, respectively. 



4. Convection driven dynamos 

Dynamos are generated by convection in rotating spherical shells for all parameter 
values as long as the magnetic Reynolds number, Rm = PmV^E is of the order 50 or 
higher and the fluid is not too turbulent where the kinetic energy density E is defined 
hy E = Ep + Et + Ep + Et- In planetary cores Pm assumes values of the order of 10^^ 
and less, but numerical simulation have achieved so far only values somewhat below 
10~^. An important feature demonstrated in figure 3 is the change in the structure of 
the magnetic field with increasing Prandtl number. While for low values of P the mean 
poloidal field is small in comparison with the fluctuating components, this situation 
reverses as P increase. This change is associated with the transition from the geostrophic 
differential rotation to the thermal wind type differential rotation caused by a latitudinal 
temperature gradient. While the magnetic energy M may exceed the total kinetic energy 
E by orders of magnitude in particular for high Prandtl numbers, ohmic dissipation 
is usually found to be roughly comparable to viscous dissipation or less in numerical 
simulations. This may be due to the limited numerically accessible parameter space, 
however. 

As long as convection does not occur in the polar regions or is sufficiently weak there, 
convection outside the tangent cylinder is nearly symmetric with respect to the equatorial 
plane as is evident from flgure 2. As a consequence the dynamo generated magnetic fields 
exhibit a dipolar or a quadrupolar character unless the non-axissymmetric components 
by far exceed the axisymmetric components of the field as sometimes happens for high 
values of Rm or low values of P. In addition dynamos of hemispherical character are 
often found, for which the field in the northern half is much weaker than in the southern 
half or vice versa. Quadrupolar and hemispherical dynamos typically oscillate. Dipolar 
dynamos often oscillate as well, but sometimes these oscillations can not be identified at 
a distance of a radius or more as is evident in the example of figure 4. 
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Figure 4. An "invisibly" oscillating dynamo in the case of P — 2,Pm ~ 1, 
T — 3 ■ 10*, 77 = 0.5, a — 1 and R — 2 ■ 10^. The plots show lines of constant in their 
left halves and meridional field lines, rsvaOdgh — const., in their right halves. The time se- 
quence of plots starts at the upper left and continues clockwise with At = 0.04 such that a full 
cycle is completed approximately. 



5. Applications to the Giant Planets 

It is obvious from the preceding sections that numerical simulations are still far removed 
from realistic descriptions of dynamo processes in the Giant Planets. There are even more 
basic issues that need to be resolved. Since the Proudman- Taylor-Theorem holds for 
barotropic fluids the strong differential rotations observed at the surface of the planets 
must be expected to continue for a considerable distance into the planetary interiors. 
Because the poloidal electric current density will be of the order UBpa where U is a 
typical zonal velocity, say lOOm/s, and Bp measures the strength of the poloidal field 
the density of Ohmic dissipation becomes of the order {UBp)^a. When this expression 
is integrated over a reasonable domain inside the planet it turns out ( Liu et al. (2006) ) 
that it exceeds easily the net luminosity of Jupiter as well as that of Saturn when the 
value of cr for metallic hydrogen is used. Accordingly Liu et al. (2006) argue that the 
dynamo is located in a region where not only the strong zonal flows penetrating from 
the surface have been truncated, but where also a is as small as is compatible with a 
magnetic Reynolds number of the order of a few 10^. Such a region would lie at a depth of 
less than 15% (35%) of the planetary radius of Jupiter (Saturn) where hydrogen is still 
a semiconductor ( [Nellis et al., 1996] ). A consequence of the latter property is a strong 
temperature dependence of a which translates into a high radial increase of A. It will be 
of interest to see how these properties and constraints can be accommodated in dynamo 
simulations producing magnetic fields similar to the observed ones. 

The dynamos operating in Uranus and Neptune are similarly constrained (Hol me fc Bloxham, 199^ . 
Traditionally an ionic electrical conductivity has been assumed for the dynamo regions 
of these planets, but recent measurements ( |Lee et al., 2006 ) indicate that water is semi- 
conducting under the relevant conditions. 

We close with some remarks on the unusual form of the magnetic fields of Uranus 
and Neptune. The properties that dipole components of the magnetic fields are relatively 
weak in comparison with other components and that the dipoles are not aligned with 
the axes of rotation are not especially unusual. Magnetic fields of such character can 
often be found in simulations of convection driven dynamos. It seems that they are more 
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Figure 5. Convection driven dynamos in the cases P = 0.3, Pm = 2, J?, = 2.5 ■ 10^ (left) and 
P = 0.5, Pm = 2.5, Ri = 2.2 • 10*^ (right). The values r = 3 ■ 10"*, ?? = 0.4 and a = 0.2 are the 
same in both cases. The plots show lines of constant Br on the surface r = ro + 0.7 (upper row) 
and r = To + 1.5 (lower row) 

often found when the diffusivity varies with radius. In figure 6 some typical examples are 
shown which seem to resemble the observed magnetic fields of the Icy Planets. 
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Discussion 

Andy Ingersoll: Jupiter's field is generated so close to the surface, i.e. at 0.86i?j. Why 
isn't the field "rougher" at the surface? 

F.H. Busse: The smoothness depends also on the effective Prandtl numbers P and Pm- 



